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A simple model to reproduce strong enhancement of dielectric response near the morphotropic phase 
boundary (MPB) is proposed. This model consists of long-range dipole-dipole interaction and composi- 
tional chemical disorder incorporated by the variation in lengths of dipole moments. By applying Monte 
Carlo simulation, we show that there appears a ferroelectric boundary phase between two types of an- 
tiferroelectric phases at an optimal strength of randomness. In the boundary phase, ferroelectric domain 
becomes remarkably large and flexible to external electric fields, leading to huge dielectric response. This 
observation indicates that huge dielectric response near the MPB originates from local polarization rota- 
tion under suppressed anisotropy by phase competition. 

KEYWORDS: relaxors, dielectric response, domain structure, morpliotropic phase boundary, Monte Carlo 
simulation 



Ferroelectric relaxors based on perovskite oxides (ABO3) 
have attracted much interest because their characteristic di- 
electric properties have several advantages to application. 
The common feature of perovskite-type relaxors is intrinsic 
randomness due to compositional disorder in the arrange- 
ment of ions on the crystallographically equivalent sites. For 
example, a typical ferroelectric relaxor Pb(Mgi/3Nb2/3)03 
(PMN)' exhibits disordered arrangement of Mg^+ and Nb^^ 
at the B sites. This intrinsic randomness is responsible for 
relaxor properties such as huge permittivity with moderate 
temperature-dependence.^ 

In material survey of ferroelectric relaxors, it is well known 
that structural phase transition induced by compositional 
change plays a special role for obtaining excellent piezoelec- 
tric properties. This phase boundary is known as the mor- 
photropic phase boundary (MPB). For example, solid solution 
(1 - x)Pb(Mgi/3Nb2/3)03-.x:PbTi03 (PMN-PT) has a rhom- 
bohedral(tetragonal) crystal structure for small(large) x, and 
shows large dielectric/piezoelectric response near MPB at 
X ~ 0.35.^ This feature is generally observed in relaxors based 
on solid solution of perovskite oxides; the MPB is located for 
(1 - x)Pb(Zni/3Nb2/3)03-A-PbTi03 (PZN-PT) at x ~ 0.09,"^ 
and for PbZri_., Ti.,03 (PZT) at x ~ 0.5.^ 

Study of the origin of large dielectric/piezoelectric re- 
sponse at the MPB is of much importance not only for un- 
derstanding of relaxor properties, but also for application to 
piezoelectric devices. Huge dielectric/piezoelectric response 
at the MPB is intuitively understood within the Landau- 
Ginzburg-Devonshire(LGD) theory.^ When different struc- 
tural phases are almost degenerate near the MPB, rotations 
of polarization become effective in response to external elec- 
tric fields instead of changing magnitude of polarization.'*'' " 
Recent discovery of a monoclinic phase near the MPB'^ is 
consistent with this polarization rotation scenario.^' " The 
LGD theory is, however, not satisfactory to describe the whole 
properties of relaxors as it treats only spatially-averaged quan- 
tities. 

In order to treat spatial dependence of polarization reflect- 
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ing intrinsic atomic-scale randomness, one of the most pow- 
erful methods is molecular dynamics (MD) simulation for an 
effective model derived from first-principles calculation. This 
approach has succeeded in quantitative evaluation of structure 
phase transition in pure ferroelectric materials, and has 
been also applied to relaxor ferroelectrics. This approach 
is, however, limited to short-time dynamics in the present 
computer resource, and is unable to access slow dynamics by 
large-scale domain motion, which is curtail to huge dielec- 
tric/piezoelectric response at MPB. In the present situation, 
statistical-mechanical approach to a simplified model is still 
valuable to clarify what the essence of physics near MPB is. 

In this paper, we propose a simple model, which can repro- 
duce several important properties near MPB.^^ By applying 
large-scale Monte Carlo simulation to this model, we calcu- 
late dielectric response, domain structure, and structure fac- 
tor Our model exhibits intrinsic competition between differ- 
ent phases, which are connected by polarization rotation. We 
show that the intrinsic atomic-scale disorder strongly modi- 
fies the shape of domain structure at MPB (transition region 
between phases). Our result indicates that in addition to the 
polarization rotation mechanism, drastic change of domain 
shape at MPB is crucial to huge dielectric/piezoelectric re- 
sponse. We also discuss its relevance to experimental results. 

We constitute a simple model to describe MPB by the fol- 
lowing strategy: (1) The model is taken so that it exhibits com- 
petition between different phases when magnitude of dipole 
moments is uniform, and then (2) atomic-scale inhomogene- 
ity is introduced by considering random distribution of mag- 
nitude of dipole moments. As a simplest model fulfilling these 
requirements, we propose a model consisting of dipole mo- 
ments on a two-dimensional square lattice: 



(1) 



Here, //, is a three-dimensional vector representing an electric 
polarization caused by the ionic displacement at site ;, and r,j 
is a displacement vector from site / to j. We further divide 
the square lattice into two sub-lattices, A and B as shown in 
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of the B-site dipole moments has spatial distribution. That is, 
we set magnitude of the B-site moment as 



A-site 
B-site 



D-AF Ferro C-AF 




Fig. 1 . (Color online) (a) A schematic picture of an arrangement of A- and 
B-site moments. Sizes of radii represent sizes of dipole moments, (b) Plot 
of energies for D-AF (dotted line), Ferro (solid line), and C-AF (broken 
line) configurations as functions of ,ub//"a. respectively. 



Fig. 1 (a), and locate dipoles with different magnitude on each 
lattice. 

Let us first consider a uniform model, in which magnitudes 
of dipole moments are constant on each sub lattice. We de- 
note the magnitudes of moments on the A- and B-site with //a 
and //B, respectively. We show the ground-state energies and 
configurations as a function of jUb/A'a in Fig- 1 (b)- assuming 
that 2x2 periodicity like in Luttinger and Tisza's preced- 
ing study. As depicted in Fig. 1 (b), there are three ground 
states. When the length of B-site moment iub/^a is close to 
unity, the ground state is a columnar-antiferroelectric (C-AF) 
state. On the other hand, the ground state is a diagonal- 
antiferroelectric (D-AF) as the ground state when //b is suf- 
ficiently small. In the jUb — > limit, the D-AF state can be 
regarded as another C-AF state on the V2 x V2 square lattice 
constructed by A-sites. It is noteworthy that directions of easy 
axis anisotropics between the C-AF and the D-AF states dif- 
fer from each other by 45-degrees. That is, the length of B-site 
dipole moment affects directions of easy axes of dipoles on A- 
site. Between the two AF states there is a wide ferroelectric 
region. The ferroelectric order has an anisotropic easy axis 
except for T - 0?^ The easy axis in this ferroelectric phase 
is vertical and horizontal near the C-AF boundary, while it is 
diagonal near the D-AF boundary. 

Analysis of the uniform dipole model has shown that the 
present model exhibits different ordered phases controlled by 
the ratio fis/f^A- We then introduce effect of atomic-scale ran- 
domness through spatial distribution of magnitude of dipole 
moments. Here, we assume that magnitude of the A-site 
dipole moments is fixed as unity (//a = 1), whereas magnitude 
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where « is an integer-valued random variable following the 
Poisson distribution P^(n) = A"e^'^/nl. The strength of ran- 
domness is controlled by setting the average value A of Pois- 
son distribution. 

Before showing analysis for the random dipole model, we 
briefly mention relevance of the present model to the MPB ob- 
served in perovskite oxides. In our model, the A-site dipoles 
corresponds to the displacement of Pb ions, whereas the B- 
site dipoles to the displacement of transition metal ions (the 
B-site in the perovskite structures). Since the intrinsic ran- 
domness is induced by compositional disorder on the B-site, 
assumption of distribution in B-site dipole moments is not far 
from reality. Then, A is related to strength of disorder, i.e., an 
compositional ratio x between a relaxor ferroelectric material 
and a pure ferroelectric material. Our model has also been mo- 
tivated from recent theoretical work focusing on importance 
of large polarization in the B-site.-^ We will discuss how the 
present effective model can be extended to more realistic one 
in the last part of this paper 

We performed Monte Carlo simulation in order to exam- 
ine effects of the B-site randomness. To realize an effective 
spin update in systems with long-range interactions, O(A^) 
Monte Carlo method was adapted. The long-range dipole 
interactions were estimated by using the Ewald summation 
method. We consider 32 x 32 square lattices throughout this 
paper To calculate dielectric susceptibilities, we apply an 
alternate external electric field with sinusoidal shape given 
as iianip sin(27rr/T) as a function of the Monte Carlo step t. 
Though we cannot define time in Monte Carlo simulation, one 
Monte Carlo step can approximately be regarded as forward 
time-evolution by a certain time step. In this paper, we set the 
magnitude and period of electric field as E^^p = 0.01 and 
T = lO'*, respectively. For thermalization, we discarded 10^ 
Monte Carlo steps while keeping on applying external alter- 
nate field. We executed 10^ Monte Carlo steps for measure- 
ment, so that 10''/t = 10^ cycles are executed for measure- 
ment. In order to average B-site randomness, 512 independent 
samples were simulated. 

In Fig. 2, we shows real part of calculated dielectric sus- 
ceptibilities per site for several values of /I as a function of 
the temperature T. As seen from the figure, the dielectric sus- 
ceptibilities has a maximum at T = Tmax, and the maximum 
temperature T^ax is lowered as A increases. The peak of the di- 
electric susceptibility is well suppressed for small A, whereas 
it rapidly grows with increasing A. The maximum value be- 
comes largest at an optimum value A ~ 5-6, and it is reduced 
for larger value of A. These features well resemble the behav- 
ior of the dielectric response near MPB in perovskite oxides, 
when we think /I ~ 5 - 6 to be the MPB in the present model. 

To examine the mechanism of emerging of the maximum, 
snapshots and structure factors S (qx, qy) for /I = 0, 5, and 8 are 
shown in the upper and lower panels of Fig. 3, respectively. 
For A - 0, the columnar-antiferroelectric (C-AF) phase is sta- 
ble at low temperatures (Fig. 3 (a)). The long-range C-AF or- 
dering is indicated by a sharp peak at {q^, qy) = (tt, 0), and 
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Fig. 3. (Color online) Snapshots of dipole configurations projected on to the lattice plane (upper panels) and structure factors (lower panels) for (a) A = 0, 
(b) A = 5, and (c) /I = 8 at T = 0. 1. Sizes of radii represent sizes of projected dipole moments, and colors of circles indicate directions of dipole moments 
on the lattice plane. Almost dipole moments reside in the lattice plane except for those located adjacent to vortices. 




Fig. 2. (Color online) Plot of real part of dielectric susceptibilities per site 
for several A's. Data of each A's are horizontally shifted to improve the 
visibility. Real temperature Tr for each A's are given by = 7-0.1(8-/1), 
where T represents a temperature given in the figure. 



most of the width of the ferroelectric domains are a unit lat- 
tice length as seen from the snapshot, explaining suppression 
of dielectric susceptibilities for small A. On the other hand, a 
relaxor ferroelectric (R-F) phase is realized at /I = 8 at low 
temperatures (Fig. 3 (c)). One broad peak of the structure fac- 
tor at (^r, ^v) - (0, 0) indicate ferroelectric short-range order- 
ing, the other peak at (q^, q^) - (n, n) reflects the diff'erence of 
magnitude of moments on the A- and B-sublattice. Here, we 
point out that at finite temperatures, direction of polarization 
in ferroelectric domain is affected by the neighboring phases. 
At /I = 8, polarization in ferroelectric domain is restricted 
along the diagonal direction affected by the D-AF phase. In 
this phase, strong randomness make peaks in the structure 
factor diffusive, i.e., extended along the diagonal direction, 
leading to butterfly-shaped structure factors which are typi- 
cally observed in relaxors. These properties characteristic of 



the R-F phase explains relatively large dielectric response for 
A> 6 m our model. 

Here, we note important difference between the C-AF 
phase and R-F phase. The former phase realized for small 
A has easy axis of dipole moments along (1,0) or (0, 1) di- 
rection, whereas the latter phase realized for large A along 
the diagonal direction (1,1) and (1,-1). This diff'erence can 
be seen by the shape of domains shown in the upper panels 
in Fig. 3 (a) and (c). We stress that transition between these 
two phases necessarily accompanies polarization rotation}^ 
What happens if A is set on the boundary between the above 
two phases? This is the main problem which we will clarify 
from now. 

The boundary between the C-AF and R-F phases is lo- 
cated around /I = 5 (Fig. 3 (b)). At /I = 5, ferroelectric or- 
dering remains in this boundary region as indicated by two 
peaks in the structure factor The ferroelectric domain size 
becomes, however, much larger as seen from the snapshot. 
The remarkable enhancement of correlation length for ferro- 
electric ordering is also indicated from the sharpness of the 
peak in the structure factor at {q^, q^) - (0, 0); the correla- 
tion length becomes much larger than the case of A - for 
which diffusive butterffy-shaped pattern is induced by local 
atomic-scale domain. A noteworthy fact is that directions of 
the polarization in ferroelectric domains are not restricted in 
contrast to the C-AF and R-F phase. This feature can be seen 
also from the shape of the domain walls, which are irregular 
shaped and much rounder than in the C-AF and R-F phase. 
The free-wheeling of the polarizing directions in the interme- 
diate region comes from competitions between the C-AF and 
the R-F states. Because of the B-site randomness, the compe- 
titions differ from area to area. In view of the local stability, 
a favorable dipole direction is affected by neighboring B-site 
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Fig. 4. (Color online) A schematic phase diagram of our model. At suf- 
ficiently high temperatures, there is a paraelectric (Para) phase. Phases at 
low temperatures depend on the average length of B-site moments /ia. Arcs 
in three squares represent favorable directions of dipoles . 



moments. However, the direction of the ferroelectric domain 
that the dipole belongs to also affects the favorable direction. 
Therefore, favorable directions of dipoles can be changed eas- 
ily by a small disturbance, so that dipoles react immediately 
to the external electric field in the intermediate phase. This 
competition is the reason why the largest susceptibility peak 
is achieved in the intermediate region. 

In Fig. 4 we summarize a schematic phase diagram of our 
model. Since an increase in size of dipole moment raises tran- 
sition temperature, phase boundary adjacent to the paraelec- 
tric phase is proportional to the average size of //b- Stable 
dipole configurations are influenced by jUs (see Fig. 1(b)), so 
that we observed the R-F, ferroelectric, and C-AF phases as 
increasing yUB (decreasing A). Due to the configurational en- 
tropy, ferroelectric phase becomes wider as increasing tem- 
perature. In the present model, the B-site randomness controls 
not only the stability of the C-AF and R-F phase, but also 
changes the easy axis of dipole moments. A schematic figure 
of the easy-axis potential on dipole moments is also shown 
in Fig. 4. In the intermediate region, the easy-axis potential is 
mixture of the ones of the C-AF and R-F phase. The result- 
ing mixed potential forms the dimple at the bottom of a wine 
bottle, and it makes dipoles easy to rotate. 

Our result is consistent with polarization rotation mecha- 
nism proposed by Fu and Cohen, which claims multiple 
ordering energetically degenerate near MPB. Their discus- 
sion is, however, based on the spatially averaged quantities 
as well as other theoretical works based on the LGD theory, 
and couldn't predict anything on domain structures. Indeed, 
our Monte Carlo simulation has demonstrated that change of 
anisotropy in dipole moments does strongly affects spatial 
structure of ferroelectric domains. This result indicates that 
local polarization rotation under suppressed anisotropy due 
to phase competition makes domain wall flexible to external 
field, leading to huge dielectric response. This is the main re- 
sult of this paper. 

Our observation of large ferroelectric domain formation 
near the MPB is expected to be examined by experiments 
such as neutron scattering^^""" and transmission electron mi- 
croscopy (TEM)."*^ ''^ Recently, Matsuura et al. observed that 
correlation length becomes large near the MPB.^^ The same 
group further fortified the mechanism by inelastic neutron 



scattering measurements and an analysis of a pseudospin- 
phonon coupled model, which is a similar scenario of our 
model. An anomalous nature of ferroelectric domain at the 
MPB has also been reported in the TEM experiment.^* 

Finally, we discuss extension of our simple model toward 
more realistic models. It is straightforward to extend our 
model to three-dimensional models without any problems ex- 
cept for computational resources. In our model, we omitted 
short-range interaction related to covalent bonding between 
neighboring ions and strain of the surrounding media. As a 
result, the C-AF phase has been realized in our model in 
the pure material side (small-/! region), instead of pure fer- 
roelectric phase usually realized in perovskite oxides. This 
artifact can be corrected by considering short-range interac- 
tion, though it is left for future work. We, however, claim 
that essence of physics at MPB should come from local po- 
larization rotation due to suppressed anisotropy by mixture 
of anisotropy potentials, which are sufficiently taken into ac- 
count in the present model. 

In summary, we proposed a simple dipole model, and ex- 
ecuted Monte Carlo simulations to it. We showed that there 
appear a boundary ferroelectric phase between the C-AF and 
R-F phases at an optimal B-site randomness, and that it has 
large ferroelectric domains with flexible walls. This novel fea- 
ture is due to mixture of two different easy axes neighboring 
phases. We expect that our model explains the reason why the 
realization of the large dielectric response around MPB re- 
gions of a wide variety of binary and ternary Pb-based solid 
solutions. 
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